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ABSTRACT 

Populations of cells that make up organ tissue grow and 
contract. A traditional approach to modeling organ size restriction 
to an observed "normal" level is to postulate a physical carrying 
capacity: effectively a limit on the physical region that can be 
occupied by the organ. The purpose of this note is to provide a very 
simple model for a cell population that grows under the control of 
positive and negative growth factors. It will be seen that such a 
model can result in logistic growth without the necessity of 
postulating a physical carrying capacity. 
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1. Introduction 

Populations of cells that make up organ tissue grow and contract in a manner 
that is roughly analogous to the fluctuations of other natural populations. Since 
organs are of bounded size their growth is not entirely uncontrolled and 
exponential overall, or otherwise a human or rat liver would either assume a 
totally outlandishly large size, or else shrink to zero. Neither such alternative is 
seen in nature, although organ sizes do vary among otherwise comparable 
members of the same species. And organs of mature hosts can change in size as a 
result of disease, toxic insult, or an operation such as partial hepatectomy, from 
which a liver can recover again to normal size and function. 
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A traditional approach to modeling population (= organ) size restriction to an 
observed "normal" level is to postulate a carrying capacity : effectively a limit on 
the physical region that can be occupied by the population. This formulation 
apparently goes back to Verhulst (1836); see Murray (1989) for recent discussion. 
In the organ situation this might correspond to a space of approximately 
pre-ordained dimension that, say, liver cells in liver tissue cannot exceed in the 
body of a mature human male. The space can be taken as given, introduced into 
other models as a parameter, and in particular cases estimated from data. It 
would be the maximum size of the liver compartment in a PB-PK model, for 
example. 

There is another alternative to the above approach that depends upon 
recognition and measurement of the presence of various biological agents called 
growth factors. There are a number of such factors that both stimulate (positively) 
and inhibit (negatively) cell population growth. Growth factors are discussed by 
Alberts et al. (1994). Aaronson (1991) provides an overview of growth factors in 
cancer; see also Rubin, Bottaro, and Aaronson (1993). The purpose of this note is 
to provide a very simple model for a cell population that grows under the control 
of positive and negative growth factors. It will be seen that such a model can 
result in logistic growth without the necessity of postulating a physical carrying 
capacity. An effective carrying capacity appears as a result of presumed growth 
factor interaction with cells. 

2. Model for a Cell Population Under Growth Factor Control 

Suppose a population of cells is of size C(t ) at time t. Its individual cell growth 
or birth rate is Xq, and its death rate is po, so its net growth rate, Ao - po, governs 
the manner and rate of growth. Starting with C(0) members, and left alone, the 
population would grow roughly like C(t) ~ which means either to 
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a large size (Ao - /jq > 0), or to zero (Ao - po < 0). Clearly such unrestricted 
behavior is inappropriate for describing a population of cells that constitutes an 
entire organ, although essentially such a model has been used to describe growth 
of tumors within an organ; see Tan (1991) for an overview; in particular the work 
of Moolgavkar and co-authors, cited in Tan (1991), is relevant. 

Now introduce a quantity a(t ) of a positive growth factor into the vicinity of the 
cell population. The amount present, a(t), changes cell birth rate to Ao + Aia(0, 
where we take Ai > 0. Also introduce a quantity pit) of negative growth factor; it 
changes cell death rate to p$ + p\ /3(f), p\ > 0. Then changing levels of a(t ) and (%t ) 
can certainly alter the properties of the cell population, from net growth to net 
decline, depending upon values of cdt) and /3(f). 

Assume that the productions of both a(t ) and /3(f) are regulated by cell 
activity in such a way that 

~ = p a C(l)S a a(l) (2.1) 

and 

PpC(t)-Spp(t). (2.2) 

That is, both are stimulated to increase by the number of cells present, and to 
diminish in proportion to their own concentration, possibly being removed from 
the cell site (organ) by blood flow or metabolism or other biological processes. Of 
course the above equations are prime candidates for replacement by others that 
more accurately depict the true interactions. 

In the presence of a and /3-factors the cells in the organ grow and decline 
according to 

= [h + *1«W]C(1) - b + (2.3) 
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So (2.1), (2.2), (2.3) form a system of three non-linear differential equations. No 
explicit solution seems immediately available, unless we make the quasi-static or 
quasi-steady-state assumption (QSSA); see Strogatz (1994) for its invocation so as 
to solve a non-linear dynamics problem along with some historical references, 
and Segel and Slemrod (1989) for a careful discussion of this approximation. 
Namely assume that a(t ) and /3(f) are able to adapt very quickly to any current 
value of C(f) to always reach a "temporary steady state": 

^ = 0 = p a C(<)- S a a(t) (2.4) 

= 0 = ppC{t) - 8pP{t) . (25) 

Adopt the approximation as true, so solve (2.4) and (2.5) for oit) and /3(f): 



a(0 = (Pa/^a) C (0 (2.6) 

and 



m = (pp/8p)c(t). (2.7) 

Let us call (p a / 8a) and (pp/ 8p) the prevalences of the a and /? factors respectively. 
Insert these into (2.3) and for convenience, put X{ = 2,] [Pa/^a)' Pi = Pl(P/?/^/l)' 

to obtain 



^ = [h + W0]c(j) - [ft, + Pic(< )]c( f ) 



= (2o-Po)C(1)-W-^')C 2 (0 (2.8) 



= (2o-Po) c (0 



1- 



Pi-Aj 

2q-Po 




This conforms exactly to the original logistic equation if the ordinary net growth 
rate, A = Xq - po, is positive, as is the effective carrying capacity 
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i/ _ ^0 AA) 



(2.9) 



Under the above conditions and starting from C(0) > 0, the population attains the 
long-run steady-state value 



C(°o) = K = 



^ Q-^0 

V\(pp/Sp)-HPa/Sa) 



( 2 . 10 ) 



The above version of carrying capacity makes intuitive sense in that 

(a) it increases with net population growth rate, Ao - po', 

(b) it decreases with increased prevalence of negative growth factor, ( pp/8p), 
and with decreased prevalence of positive growth factor, (p a / 5^); 

(c) the inhibition effect of negative growth factor, p\ = p\ (pp/5p), must 
exceed the stimulative effect of the positive growth factor, 

A-i = A,] {pa/ <5a)- 

If any of the above conditions are violated the population development 
becomes radically different, but can also be interesting. 

The time-dependent population size is seen to be of the familiar logistic 
growth form 



KC(0)e 

~ X-C(0) + C(0)c Af 



with K as in (2.10), A = Ao - po > 0, and 0 < C(0) < K. 

Note that the formula has biological meaning even if C(0) > K, and also if 
K < 0: suppose that A = Aq - po > 0 but X{> p{; then write K' = -K > 0 to get 



= ac(j)[i+c(i)/k'], 

the solution to which is 



( 2 . 12 ) 
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(2.13) 



C(t) = 



jC'C(0)e 



At 



K'+C(0)-C{0)e 



At 



if t< — In 
A 






1 + 






K' 

C(0) 



1 



; it explodes when f = — In 



1 + 



K' 

C(0). 



; this might plausibly 



model an especially malignant tumor growth. Finally when A = Ao - po < 0 and 
K' = -K > 0 we simply get (2.11) with -A replacing A, once again a logistic model, 
but now one that decreases as t increases. 



3. Stochastic Models 

The above model can be "made stochastic" in various ways. One is to re-state 
the growth factor and cell-growth equations as a system of three non-linear 
Ito-type stochastic differential equations. Analytical solutions are not likely to be 
available, but some asymptotics might well produce explicit results. 

Another approach is to assume that a stochastic version of C(f), namely C(f), is 
a birth-and-death process with transition rates copied from the right-hand side of 
(2.8). For example, C(t ) evolves over state space (0, 1, 2, ...) according to 

p{C(f + dt) = C{t) + l|C(f) = /} = (Aq - fi 0 )idt + o{dt) (3.1,a) 

P{C{t + dt) = C{t)-l\C{t) = i} = (ii[-Xi)i 2 dt+o{dty > (3.1,b) . 

these holding when A = Aq - po > 0, 0; otherwise modification is needed. 

Although E[C(f)|C(0) = C(0)] * C(t) because of the non-linearity in the generator, 

(3.1), it can be of interest to study the above stochastic version's transient 
properties, such as first-passage times from low states (population size) to high, 
or the reverse, e.g. to C(f # ) = 0 when the population dies out. 

The described approach essentially minimizes attention to the stochastics of the 
growth factors and ignores non-linearity, hence is a prime candidate for an 
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upgraded treatment. Nevertheless it is appealing for its simplicity and easy 
availability, and is offered as an interim approach. 

4. Summary 

It is shown that classical logistic growth can be induced in a non-traditional 
manner by hypothesized action of growth factors, rather than by action of a 
physical carrying capacity (although the latter may operate as well). Modification 
of the effective carrying capacity to be negative has biological interpretability. 
The resulting models may perhaps find a use in cell proliferation and cancer 
modeling. 
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